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ON THE FORMATION OF THE DERIVATIVES OF THE LUNAR 
COORDINATES WITH RESPECT TO THE ELEMENTS* 

BY 

ERNEST WILLIAM BROWN 

In the determination of the action of the planets on the moon, it is usual to 
use the method of the variation of arbitrary constants. This method requires 
the knowledge of the derivatives of the lunar coordinates with respect to the 
angular elements e, it, the epochs of mean longitude, perigee and node — and 
also with respect to the other three elements n, e, 7 — the mean motion, eccen- 
tricity and sine of the orbital inclination. If we have a theory in which the 
coordinates, so far as the actions of the sun and moon only are concerned, are 
expressed literally in terms of these six elements, it is a simple matter to find 
the derivatives. 

But this is practically not the case. The convergence of many of the coeffi- 
cients when expressed in powers of m = n/n' — the ratio of the mean motions — 
is so slow that a literal theory with the required accuracy seems almost impossible 
on account of the labor required for its development ; the slowness of converg- 
ence does not occur as far as powers of e , 7 are concerned. The theory of Hansen 
is entirely numerical and therefore only the derivatives with respect to the angu- 
lar elements can be obtained from it. The theory f which I hope shortly to 
bring to a conclusion is semi-numerical, i. e., the numerical value of m is inserted, 
the other constants being left in a literal form. J 

In considering how this theory was to be adapted so as to be of use when con- 
sidering planetary action, the only difficulty was, therefore, the formation of the 
derivatives with respect to n. It is this difficulty which has been solved in the 
following pages. The semi-numerical theory once completed, it is shown that the 
derivatives with respect to n can be made to depend on quadratures with respect 
to the time. The formulae are put into a form ready for computation and it is 
shown how a simple transformation will permit us to use the advantages of can- 
onical sets of constants without being hampered by their disadvantages for pur- 
poses of calculation. 

* Presented to the Society February 28, 1903. Received for publication April 14, 1903. 

t Theory of the Motion of the Moon, Memoirs of the Royal Astronomical Society, parts 
I, II (1897, '99) in vol. 53 ; part III (1900) in vol. 54. This will be referred to below as T. M. 
M., parts I, II, III. 

t A detailed statement will be found in arts. 4-6. 
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§(i) Transformation of the integrals of the equations of variations. 

1. Let cc,2/,a,x,^, 2;,be the coordinates and velocities of a particle moving 
under the action of a force function F ., so that 

(1) 

the axes being fixed. Denote by a (p = 1, 2, • • •, 6 ) the arbitrary constants 
of this motion and put 



d'x dF 


d^y dF 


dh dF 


dt^ ~ dx ' 


df - dy ' 


df~ dz 



da ' p dtda da dt da ^ ' ' 

i> p p p 

where the solution is supposed to have been found and expressed in terms of the 
arbitrary constants and the time. It is more convenient at the outset to take 

the a to be constants and not some of them to be variables of the form b t -j- B 

p p p 

as is usually done in the applications ; certain changes to the latter form will 
be made later. 

The equations of variations are 

d'-F d'-F d"F 

^^> "'^"" 5x^*^ + 5^5^^"+ axa^^^' 2/p=---' %=■■■• 

They possess the 15 integi-als, 

(3) (p, 9) = x^x^ - xx^+ y^y^-y^Vp + %^^ - ^^% = Gj,^ (p, 9=1, 2, ■ • ■, e), 

where C is a constant and C = — C , C = . 

PI PI IP pp 

Let the solution be such that x, y , z each consists of sines or cosines of sums 
of multiples of the angles 

^,-<+«- (i = 4,5,6), 

where h. and the coefficients depend only on the arbitraries ctj, a^., a^. The 
solution may also contain other known angles and constants present m F; it is 
unnecessary to specify these here. Then it is well known that 

(4)C,,= 0, C^,, = 0, C,.= '^, 6,.= --|? (M'=l,2,3;i,/=4,5,6), 

• i-3 

where c,, c^, Cj are functions of ctj, a^, a^ and the known constants, and B is a 
constant expressed in terms of c^ , c^ , c^ (instead of a, , a^, a^ ) and the known 
constants. 

2. Consider the determinant of the sixth order 

^ = \x^,x^, y^,y^, z^,z^\ (p = i,2, •••,6), 

and form dAfdt. The result contains six determinants obtained by differ- 
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entiating the elements in each of the six columns. The determinants formed 
by differentiating the elements in the second, fourth and sixth columns will van- 
ish owing to their possession of two identical columns. Those formed by differ- 
entiating the elements in the first, third and fifth columns will also vanish when 
we substitute for x , y , z their values from equations (2). Hence dA/dt = 
and 

A = constant = JT. 

To find IT write A in the form 

and multiply the two forms for A by the ordinary rule for the multiplication of 
determinants. Using the equations (3) it is evident that we obtain 
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a comparison of terms showing that the positive square root is to be taken. 

Let k.. denote the first minor of C.. in this determinant with the sign so taken 
that a cyclical interchange of the first three or of the second three suffixes does 
not require a change of sign in the minor. Then 

(5) 2.^|..C^,. = A or according as i' = i or 4= *• 

3. Next, denote by A the first minors of the elements in the first column of 
A and by A the second minors of those in the first two columns. Then 

J pq 

(;, = 2,3, ••■,6) 

by (3) 

= ^i i %^ Vp^ y^^K^%\- ^42^42 - G,^A,^ by (4). 

The first determinant in the right hand member is, as before, easily seen to be 
equal to d\jdt. Hence, writing down the two similar equations, we have 



x^A, 



I 4 p 



Vp ^Vpi ^p-^'^pS 



x^\ 

x.A, 



«4^1 = - ^^42 ^42 



^43^43' 



■ «. A, = — C„A., — C-,A,,, 



x,.A, — x^A, 



^62 62 ' 



^63 63" 
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Multiply these equations by k^^, h^^, k^^, respectively and add. The right 
hand member vanishes by (5) and the left hand member becomes integrable, 
giving, since we have similar forms for A^, Aj, 

A. = (k..x. + k.,x. + k-.x^^A. (t-^l, 2, 3), 

where ^^ is a constant. Similarly 

A^, = ( \. x^ + \. x^ + ^^3^. X3 ) A. (i = 4 , r, , 6 ) . 



Put 
Then 



X = k,^x^ + k.^x^ + k^x^, Xj = \.x, + \.x^ + \.x^. 

l\=^.A.X.X.— ^.A.X.X. (» = l,2,3;i=4,5, 6). 

To determine the constants A , we notice that if any two suffixes be inter- 
changed in the determinant A , it merely changes sign ; hence A.= A.^== 1, 
and 

A = ^.x.X.-^.x.X., 
which may also be written 

^.{x.X. — x.X.)^K (f=l,2,3). 

Again, by the properties of determinants, 

= Scc A =2v A ==1.11 A = 2s A = 2 2; A (p=l,2, ••■,6), 

P p Op p up p P P P P III 

so that, if 1^., Z^ be defined with reference to t/., z. in the same way as X ^ was 
with reference to «., and if we write down the corresponding expressions for A 
obtained by interchanging x, ?/, a, we have the following set of 15 equations: 

^ = ^{x.Y.-y.X.) = 1(x.Z.-z.X.) 
= -2(x.F.-y.X.) = I.(x.Z.-z.X.), 

with the eight other equations formed by a cyclical interchange between x,i/, z; 
the summations refer in all cases to the values i = 1, 2, 3. Of these fifteen 
equations, not more than ten are easily seen to be independent of one another, 
whatever be the values given to x^, y^, z., X^, Y., Z^. 

§ (ii) Statement of the problem and its form,al solution. 

4. The main object of this paper is to determine the derivatives of the coor- 
dinates of the moon with respect to the mean motion m, it being supposed that 
the derivatives with respect to the other arbitraries have been obtained. The 
results found above are true whenever the coordinates satisfy the specified con- 
ditions and the same may be said of most of those which follow ; but it will be 
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convenient to put them into a form which will be applicable directly to the 
lunar theory, since they are required for the determination of the action of the 
planets on the moon 

It is well known that the coordinates of the moon referred to fixed axes are 
expressible by means of the five angles 

(8) ?i< + e, n't + e', cnt + e — tt, gnt -\- e — 6, n't-j-e' — tt', 

where n , n' are the mean motions of the sun and moon ; (1 — c)n, (1 — g)'>^^ 
those of the perigee and node of the moon; e, tt, 6, the usual angular con- 
stants of integration referring to the moon ; e', tt', the known angular constants 
of the sun's motion. Also, c , g and the coefficients are expansible in powers of 

n , a 

]^'=^m, e, 7, e, -^7 = «, 

where e, 7 are "eccentricity" and "inclination" constants of integration refer- 
ring to the moon ; e' (eccentricity), a' (mean distance), known constants of the sun's 
motion ; n^a^ = sum of the masses of the earth and moon ; n' a' = sum of the 
masses of the earth, moon and sun. The set of arbitrary constants a^, a^, a^ is 
taken to be w , e , 7 and the set a^, a^, a^to he €,-7r, 0; n,n[l — c),«.(l — g^ 
are the motions of the corresponding angles. 

If the coordinates be referred to axes of which one points to the mean place 
of the sun or moon, the coordinates are functions of four angles only, the first 
two of (8) being replaced by their difference. 

5. When it becomes a question of finding the derivatives of the coijrdinates 
with respect to the arbitraries, there is no difficulty as far as e, tt, ^ are con- 
cerned, since these only occur under the signs, sine and cosines, whether the 
theory has been worked out literally or with numerical values for m, e, 7, e', a. 
But this is not the case with », , e , 7 : a literal theory like Delatjnay's is nec- 
essary. But the use of such a theory gives rise to a practical difficulty. Every 
coefficient of a periodic term is of the form aA\ where X is the product of 
the highest positive powers of e, 7, e', a which are factors of the coefficient, 
and A is expansible in powers oi m, e^ , 'f, e' , a^ . If we arrange A accord- 
ing to powers of the last four parameters, the convergence is sufficiently rapid 
for practical purposes and if, further, the coefficient (a function of m only) 
of each of these powers is determined with sufficient accuracy, there is no diffi- 
culty in finding the derivatives with respect to e , 7 . It is quite different when 
we arrange according to powers of »i . In the first place, the rate of conver- 
gence is in many cases very slow, corresponding frequently to that of a geometric 
series with a common ratio J, and the rate of convergence of the earlier terms 
is less in the series formed by taking the derivatives with respect to n or 
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m . The most accurate literal theory is that of Delaunay and this is insuf- 
ficient for modern demands. Moreover, it may be said to be almost a practical 
impossibility to calculate a literal theory to a sufficient degree of approximation. 

6. It was with these facts in view that I set out to obtain a semi-numerical 
theory in which the parameters e, e' , 7, a were left in their literal form, while 
the numerical value of m (which is known from observation with very high 
accuracy) was used from the start. The values of the coordinates have been 
published * as far as the fourth order inclusive in these four parameters and 
this order is probably sufficiently high to obtain the derivatives of the coordi- 
nates with respect to e, 7 with the degree of accuracy required to obtain the 
action of the planets on the moon. But we cannot, of course, directly obtain 
from it the derivatives with respect to m. It might be possible to use Delau- 
nat's theory in conjunction with the numerical results by a somewhat empirical 
method which I adopted in an earlier paper f and which served well the purpose 
for which it was intended. But it appears doubtful whether it would be suffi- 
ciently accurate for the determination of all the planetary inequalities. I have 
therefore considered this problem : 

Given the derivatives of the various functions with respect to e, tt, 6, e, 7, 
to find those with respect to n from a theory in which the numerical value of 
m has been substituted. 

The possibility of solving this problem practically arose from considering a 
subsidiary result which I gave in an earlier paper, J where it was solved so far 
as the terms independent of e, e', 7, a were concerned; the object and the 
method were, however, quite different from those of this paper. In the follow- 
ing sections the problem will first be solved in the form which naturally presents 
itself and afterwards the results will be so transformed that certain difficulties 
arising from the infinities of the integrals are avoided and that the results of 
the theory, above referred to, can be directly used. 

7. It has just been seen that we require to determine one set of the deriva- 
tives in terms of the other five : let this set be x^ , y^, z^. For this purpose we 
have the fifteen equations (6), (7); we have seen, however, that they are not all 
independent. Let us attempt to determine X-^ without an integration from 

l{x,X,-x,X,) = Jr, ^{x,r,-y,X,)=0, ^{x,r,-y,X)=^0, 
of which the terms containing the suffix 1 are 

*l-^l '*'|-^i? "^l-"-! 2/1-^1' ■^'l-'l 2/1-^1' 

*T. M. M., parts I, II, III. 

t On the theoretical values of the secular accelerations in the lunar theory, Monthly Notices of 
the Royal Astronomical Society, vol. 57 (1897), p. 346. 

i On the solution of a pair of simultaneous differential equations, etc., Memoirs of the Cam- 
bridge Philosophical Society, vol. 18 (Stokes Memorial, 1900), pp. 94-106. 
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and the capital letters contain the suffices 4, 5, 6 only. The Jacobian of these 
with respect to x^, x^, y^ is zero ; a similar result appears to hold if we attempt 
to combine any of the equations. Hence we shall use the first two of (6) to 
determine aSj , y^, and one of the others to determine z^ . We find 

(9) x, = X,f ^- -.^^^-.^^^ ±-3^^,, 

with a similar equation for y^ ; no added constant is necessary since x^ cannot 
contain terms of the form AX., that is, of the form Ax., for the x. consist of 
sines and x^ of cosines. It will be sufficient to consider the equation for x^ 
since that for y^ is exactly similar ; the equation for »j is treated separately. 

8. There are three points to be considered in connection with the equation 
(9). First, ajj, x^, x.^, x^ give rise to terms containing the time as a factor owing 
to the presence of e , 7 in c , 17 , that is, in the arguments of the periodic terms ; 
practical convenience requires this to be avoided. Second, JS^, X^ contain the 
derivatives dc./dn owing to their presence in certain of the coefficients k^. ; these 
must either be calculated or eliminated. Third, the principal term of X^ , 
which appears as a denominator under the integral sign, is the derivative of 
a cos (wi + e) with respect to n and it may therefore vanish. In order to show 
how the first two difficulties may be avoided I shall develop the results with a 
canonical system of constants and then adapt them for practical use by the adop- 
tion of a certain semi-canonical system. The third difficulty easily disappears 
when complex coordinates are used. 

9. Instead oi a^ = n, a^ = e, oij = 7 as the first set of arbitrary constants, let 
us take ctj = Cj , a^= c^, 0:3=03, and suppose that it be required to find the 
derivatives with respect to Cj, those with respect to c^, c^, e, tt, 6 being sup- 
posed known. We have 

J^=i, hi+s = 1 ' ^ij = ( J 4= j + 3 ) , 



X. = x.,,, Y. = y..„, Z.=. 



'i+3> 



and the first of equations (6) becomes 

(10) 2;.(x.x.+,-fl;.i.+3) = l (i = l,2, 3). 

Denote the three angles by w.= b.t + /8. (;' = 4 , 5 , 6 ) . Then 

dx dx dh. 

dc. dc. ' dc. -i' 



3 d dx / ^ \ ( •<? '^^j \ ^ dx '^ ( -r (^'^j \ 

:. = ^« a^ + V *dt)\ Jdc.'^Jj^d^. di '^ *dt\ ''d~c.'^')' 



dx d dx 
dc 
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where ^/^c^- denotes that the derivative is taken with respect to c^, only in so 
far as c. occurs in the coefficients of the periodic terras of the function on which 
it operates. 

Substitute in (10) and equate separately to zero the terms which do and do 
not contain i as a factor. The latter give the well known results 



The former give 



or 



< 





dbj 
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db. 
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) = 1, 




, - x.^^x. 
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^.-2, 
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(j=l,2,3;J = 4,5,6). 



^t+igf'^i *.+s'*'.- I + -^/^i J„ ^i^i+x — •'■' 



according as we use the second or first form of dx/dc.. It is necessary to use 
the first form for i = 1 on account of the integration with respect to t and it is 
better to use the second form for i = 2 , 3 . The equation then becomes 

d dx . dx dx . dx dx . dx d.h, 

^*dtdc^ ""'^dc^-^ ""'dc^^-^^'dc^ ''«5c3+*''ac, "^^ Jdc,^- 

, db. 



= Q- 






dc, ' 



where the definition of Q is evident. 

Dividing by x\, putting db^jdc^ = dbjdc.,, db^/dc^ = dbjdc^, and integrat- 
ing, we obtain 

The only derivative with respect to Cj in the right hand member is db^/dc^ . 
But since dx/dc-^ can not contain any term factored by t, the constant term 
under the integral sign must vanish. This fact gives 

db, . ■ ,. Q 

-y- = const, term m the expansion of— \, 

1 4 

and the determination of dxjdc^ is complete. 

10. The objection to leaving the results in this form for practical application 
is the fact that the formation of the coordinates as functions of c,, c^, c^ is very 
troublesome, especially as far as c^ is concerned. But I shall show later that 
there is no objection to the use of c^ , c^ instead of 6,7. Hence it will be con- 
venient to develop the results for the semi-canonical system J^ ( = w ) , c, , Cj , so 
that Cj is considered to be a function of the independent constants w , c^ , c^ . 
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With this semi-canonical system, we find 

^ dc^ „ dc^ ^ dc.^ db^ dc^ db^ dc^ 

'■* dn ^* dn '* dn dn dc^ dn de^ 

^=^'25 = Cjg = ^14 , = C,5 = Cgj = Cjg = C,j = ^24 = he — hi = K5 ' 

J J -j-r' tie, -I tiiC', J etc. 

■y^ LvL', etc. "YT" WC, T^ etc. 



dn dn 



, ^/ , t?6, db, \ 5x , ^„ 



C&iC Oj O ''JO XT' , Vt" 

dn dt dn ^ ''^ 

and the new form of equation (9) is 

where X^, has the value just given and 

There are no derivatives with respect to n in (11) except dc^ld?i and this is 
determined by the vanishing of the constant term under the integral sign. It is 
remarkable that the expression for Q' which arises with the semi-canonical system 
is shorter than that for Q which arises with the canonical system. The equa- 
tion (11) might have been obtained from (10) by direct transformation of the 
set Cj , Cj , Cj to the set n, Cj , Cj . 

11. There are two determinations to be made. The first is to find the values 
of c, , c„ , Cj from a theory expressed in terms oi n, e, y and the second is to ex- 
press this theory as well as c^ in terms of n^, c^, c^. 

For the first, I have shown * that 

C; = const, term in ajx,^, + yy.^^ + zz.^^, 
which permits of the calculation of c. with ease. We find 

c, = ■^/aMA(e^ 7^), c^ = eWafiI'^{e',y'), c, = 7' 1/ a/tP, (e% 7' ) , 

* On certain Properties of the Wean Motions, etc., Proceedings of the London Mathe- 
matical Society, vol. 28 (1896), p. 150. The proof involves the supposition a/a' = 0. That 
the result holds when aja' + will be shown in a paper to appear later in the Transactions. 
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where P^, P^, P^ are positive power series in e^, 7^, and jj. is the snm of the 
masses of the earth and moon. It is understood that the numerical values of the 
other constants may have been substituted ; it is unnecessary in any case to 
specify them. From the second and third of these equations, the values of e^, 
7^ are quickly deduced as power series in cJV' afi, c^li/a/i, particularly as it 
will probably never be necessary to go beyond the fourth powers of e , 7 . Hence 
c^ is easily expressible in powers of the same quantities and likewise 6j.. b^. 

The coefficients contain odd powers of e, 7 in general and therefore will in- 
volve powers of the square roots of the new parameters. The final transforma- 
tion is therefore to put 

C2 = E^l/a/ci = E^a^n, c^= r^\/ afji,= T^a-n, 

so that in Q or Q' we can put 

Bid did 

dc^~ 2a^nE Be^ Bc^~ 2a^nr Bt^ 

where the coordinates are transformed by the equations 

e = EP.^{E\r^), v = rP^{E\T''). 

It is to be noticed that the value of z^ may be obtained without the necessity 
of an integration from either of the equations 

(12) ^.(x.Z. - z.X.) = = 2.(y. Z - z. Y.), 

\ / l\tt t X J %\i^ % t 11/' 

when x^ , y^ have been obtained. 

12. We have next to transform to complex coordinates to avoid the occur- 
rence of vanishing denominators. Put 

u = x + yV—i, u. = x. + y.-i/—l, U^= X.+ r.V— i, 

s = x — yV—l, s^ = x^ — y^^/—l, /S'. = X. — Z: l/— 1 . 

from the equations (6), (7) we easily obtain 

(13) = ^J^u,U,-u,U,) = ^,{u,Z,-z,U,), 

(14) ^^{u,S,-s,U,)='iK, ^^{n,S,-s,U,) = ^,iu^Z^-z,UJ = 0, 

with the corresponding equations formed by changing the sign of y — 1 . The 
first of (13) determines w, , with the denominator U^ which does not vanish * ; 
the second equation then gives », . The other equations may be used for veri- 
fication purposes. 



* Since X^, T^ never vanish simultaneonsly. 
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It is quite evident that, in order to obtain the equations corresponding to 
those of arts. 7—11, all that is necessary is to substitute m for x, Uiov X, and 
omit the term unity in Q and Q' . To find dzjdu we have 

c5s Mj 1 dc^^du dz du dz \ 

^ ^ ~&u^V[ '~Tr^dn\ dc^ ^' ~ 53, "'^ "^ d^^ ''<' ~ d^^ '"V ' 

where 



dz Uy ry \ dc^ ( du 
du~ U, ' U,~dn\dc.^-'' 


dz du dz 
5c/''^ + ac, -'^^ dc^'^^ 


dc^ dc^ 


dc^ cZc, 

^^-'''~dc/'^~dc^ 



3 



13. As a rule it is more convenient to have the coordinates referred to axes 
of which one is directed to the mean place of the sun or moon. There is no 
trouble in either case in making fhe change. We substitute 

M exp. + i/— 1 (?i'^ + e') or w exp. + i/— 1 («^ + e), 

in the respective cases, for u. The various small charges of this nature which 
experience has shown to be best for actual calculation are made in the summary 
which follows. It is scarcely useful to give all the details of the algebraical 
work which leads up to the results. The formulas are somewhat long and neces- 
sarily unsymmetrical, as is generally the case when they are best adapted for 
numerical computation. It is to be noted also that the calculations are of the 
same nature as those which constitute the bulk of the work in my lunar theory, 
namely, multiplications of series. The great majority of these multiplications 
are, however, already done, having been previously needed, so that the new work 
is much less extensive than would appear from the formulae. The order of accu- 
racy with respect to e, 7 to which they give the values of the derivatives with 
respect to n is one less than that of the coordinates, but in any case it can be 
found to a sufficient degree of accuracy for practical purposes. 

§ (iii) Summary of the results adapted for calculation. 

14. The notation adopted in this section is that of my theory.* It is there- 
fore necessary to substitute for the symbols denoted above by u, s, the values 

u exp. V' — l(n't + e'), s exp. — V — l{n't + e' ) 

to reduce to axes moving with the sun's mean motion. We shall suppose that 
this substitution has been made. 

Now u is expansible in powers of f, f", ^^, ^°'. Let 

d d d d 
n = f ' D — t" T) —t^ T) — ?■" 



*T. M. M., parti. 



1903] OF THE LUNAR COORDINATES 245 

a notation which can be adopted since c, g, ni, are incommensurable with unity 
and their numerical values are not substituted in the index of f . Then 

Put* 

a.u' = u^~'^ . as' = st, as' = s; 

then w'a is the complex coordinate referred to the mean place of the moon as 
axis of X . 

15. I shall now indicate briefly the various steps required to find the deriva- 
tives with respect to six arbitraries. 

( a ) The derivatives with respect to e , tt , ^ are given by 

du ^ , , -r dii 



( 6 ) The canonical constants Cj , c^ , c^ are given by the method referred to in 
art. 11. The equations for them and for certain new constants (not quite the 
same as those of art. 11) introduced for convenience, are 

(n - n )a-A = Cj = - 1(m - n )a2 [(Z> + 1 + m)^*' • (Z>j + D^ + Z>^ - 1 )s' 

+ (Z> - 1 - m)s' ■ (Z>, + D„ + Z»^ + 1)m' 

(_B) _(n-n')a2E^ = C2 = -|-(w-?i')a^[(Z>+ 1 + m)^' • D^s' 

+ (Z»-l-m)s' ■ Du' + 2Dz' D/y , 

— (w — 7i')a"r^ = Cj = |(w — w')a" [(Z> + 1 + m)-i(' D s' 

+ (Z>_l_m)s' -Bu'^-IDz ■ B^z'Y, 

where [ ]° denotes the constant term in the development of the enclosed func- 
tion and /x is put equal to unity. The actual calculation from these rather long 
formulae is really quite brief. The operators Z)^, D^, D introduce only in- 
tegers and the necessary multiplications of series are nearly all at hand; with 
these, a few hours' work will suffice to determine Cj, c^, Cj to the fourth order.f 
We thus find A , E^ , r^ expanded in powers of e^ , k' (the constants of my 

*The symbol u' is not used in T. M. M. 

fThe relations of a to Del aon ay's final L, G, H, are given by 

c, = L, c, = {G~L), c,= {H~G), 

wheni«=l. 

Trans. Am. Math. Soc. 17 
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theory) with numerical coefficients if the numerical values of e' , a have been 
substituted, and in powers of e' , a^ if their numerical values have not been sub- 
stituted. 

( c ) Revert the series for e" , r^ so as to find e^ , k^ expressed in powers of E^ , 
r^. Then find A, c, g, in terms of the latter constants. Take the positive 
square roots to obtain e, k. 

((?) Substitute the values of e, k in the expressions for u , s' , so as to obtain 
the coordinates in terms of E, r. 

(e) Put 

^ E 5e e 5e k 5r^ e K dv^ 

dc, 



^ L2Z7'^J' 



dn = -^^^ 



and compute V , Q and the constant /3 . Then we find the derivative with 
respect to n from the formula 

no added constant being necessary on performing the operation D~^ . The 
operations dJdE, d/dr are performed only on the coefficients of the various 
powers of f in u' , Du' . The final result is the value of d(u' n) jdn when E, r 
have been replaced by their values in terms of Cj, c.^: the change was a tem- 
porary one for convenience only. 

The actual calculation of ^/C/'"is performed in the following way. The 
principal terms in u (those depending only on m and therefore of the form 
S.a. f^' where the a. are numerical coefficients and «(,= !) being denoted by w^ , 
the principal terms in U' are given by 

Hence 

1 1 f^ ^u'-m ^(U'-u' 



U''~ w' 



i-^^;^-r^y- 



The functions 1 / U'^ , being power series in f '^ are best calculated by the 
method of special values. It is to be noted that all terms on which D^ has 
operated contain the factor e , and those on which D has operated, the factor r^ . 

In the expansion oilj U' , the value of u' is used as far as the order actu- 
ally required; in that of Q, the value of w' to one order higher in E, r. 
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The values of du' /ds, du' jdv have been obtained In the course of the inves- 
tigation and we therefore have the derivatives of with respect to all six con- 
stants. 

(_/■) Finally separate into real and imaginary parts and obtain the derivatives 
of the two coordinates referred to moving axes, that of x pointing to the mean 
place of the moon. For the transformation we have 



t' = cos iD -f V— 1 sin iD , ?'" = cos il -f l/— 1 sin il , 

(D) 

^'^ = cos iF -f l/— 1 sin iF, t™ = cos il' -f j/— 1 sin il\ 

where D, I, V , F are the usual angles in Delaunay's notation. 
{g) Put 

Z' = (Z», -M) + Z> )s' - ^„ i» 2i' - ^2 Dz\ 
and determine dzjdn from 
(F\ n-n'dz ( /3Q \ B 

where 

^, du' D z dz D u' ^ du ^ , dz D„u' 
^ ~ De e 5b e "^ ar^ ^K- ar r ' 



The expansion of If W is made as before and the calculation of the various 
functions follows a quite similar plan. 

16. The value of d(^u'^a.) /dn has already been obtained. If the suffix zero 
denotes the parts depending only on m , the connection between the results here 
and those in T. M. M., pt. II, chap. V, section (ii), is given by 



I have found elsewhere * by a totally different method, 
which gives the same value for c?(Cj)j/(?n. 



*Page 153 of the paper referred to in art. 11 and page 348 of the first paper referred to in 
art. 6, 
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18. If it be intended actually to use the canonical equations for the variation 
of arbitrary constants, namely, 

dc. d R dw. dR 

di=a^; 'dt=-B^+^^ (i=l,2,3;j = 4,5,6), 

it is easy to make the necessary changes. Denote by parentheses the derivatives 
when Cj , c^ , Cj are considered to be the arbitraries. Then if (^ be any one of 
the functions used, the derivatives of ^ with regard to c^, c^, Cj are found by 
solving the equations 



\dc^^)dn' 'd^^~\dcjdc^'^\dc^,)' 'd^" \ dVJ d^^'^ \ dcj ' 



8(f> / d <p\ dc^ d <p 

dn 



with respect to (5^/5cj), {d^jdc^), {d<fildc^y, dcjdn= —/3a'' is given by 

art. 15 (e) and dc^/dc^, dc^/dc^ will have been found. 

Haveefoed College, 
April 14, 1903. 



